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We study the SU(3) gluon propagator in renormalizable R^ gauges implemented on a symmetric 
lattice with a total volume of (3.25 fm)^ for values of the guage hxing parameter up to ^ = 0.5. As 
expected, the longitudinal gluon dressing function stays constant at its tree-level value Similar 
to the Landau gauge, the transverse R^ gauge gluon propagator saturates at a non-vanishing value 
in the deep infrared for all values of ^ studied. We compare with very recent continuum studies and 
perform a simple analysis of the fonnd saturation with a dynamically generated effective gluon mass. 


Introduction. Dynamical mass generation is arguably 
one of the most important non perturbative features of 
QCD. Within the quark sector, this phenomenon origi¬ 
nates from a cloud of low-momentum gluons attaching 
themselves to the current quark. At the theoretical level 
this effect is described by a so-called (quark) gap equa¬ 
tion; its solutions describe the evolution of a (chiral) cur¬ 
rent quark of perturbative QCD into a constituent quark 
(with a mass ^ 350 MeV) as the momentum flows from 
UV to IR values. This dynamically generated mass ac¬ 
counts alone for the largest part of the proton’s mass, 
explaining how hadron masses emerge (dynamically) in 
a universe with light quarks. 

In the gauge sector, the propagation of the aforemen¬ 
tioned gluons is also described by a gap equation. In the 
Landau gauge, its solutions saturate in the IR to a non¬ 
vanishing value. This can be interpreted, as in the quark 
case, through the concept of a dynamically generated 
(and momentum dependent) mass [TH5]. The associated 
gluon mass function m? is a monotonically decreasing 
function that appears in the transverse part of the prop¬ 
agator (hence gauge invariance is not tampered with); 
in addition, it is power-law suppressed in the UV (and 
therefore unobservable in perturbative applications). 

The dynamical generation of a gluon mass substan¬ 
tiate the hope that QCD is nonperturbatively well de¬ 
fined, implying that the theory generates by itself the 
needed IR cutoffs. In addition, its presence affects the 
theory in a variety of ways: for example, it gives rise to a 
quark-gluon interaction coinciding with that required to 
describe ground-state observables m- 

Lattice gauge theories have been instrumental in un¬ 
equivocally establishing IR saturation in the gauge sec¬ 
tor ElEHm. In fact, once discretization artefacts are 
accounted for, their results are equivalent to exact all¬ 
order results, and therefore are regarded as benchmark 
tests for the different continuum approaches OIIBHSH]- 
While, in general, lattice computations do not need any 
gauge fixing (GF), the latter becomes essential for the 
computation of Green’s functions (propagators and ver¬ 
tices), as it is well-known that these objects depend on 
the gauge condition employed. This is not a problem per 


se, as there are many instances in which physics is bet¬ 
ter understood in a particular gauge {e.g., partons make 
sense only in the light cone gauge); however, it clearly 
highlights the importance of performing simulations in 
as many gauges as possible, in order to discern which as¬ 
pects of the nonperturbative behavior of the function at 
hand are (or are not) affected by a gauge choice. 

While the gluon two-point function has been studied in 
covariant and non-covariant gauges |39fl42| . reliable cal¬ 
culations in renormalizable-^ (R^) gauges |13] have not 
been systematically pursued so far. Despite being the 
only class of gauges which are completely under control 
at the perturbative level, its lattice implementation has 
in fact proven to be quite complicated |44lH5] . due to 
a no go result |44| . A viable lattice formulation was fi¬ 
nally put forward in Ref. m- In practice, however, 
one still encounters significant convergence problems m]- 
I53| when attempting numerical GF, which unfortunately 
become more severe as the GF parameter ^ and/or the 
lattice volume become larger, and the number of colors 
Nc and/or the lattice coupling /3 are small. As a result, 
there have been only preliminary studies of the gluon 
propagator 

In this letter, we present the SU(3) gluon propagator 
in gauges for a relatively large lattice volume (3.25 
fm)'* and a GF parameter up to ^ = 0.5. This allows 
us to address in some detail the IR behavior of the 
gluon propagator, and study the dynamical generated 
gluon mass beyond the Landau gauge limit. 

R^ gauges framework. In the continuum, gauge fixing 
is achieved by adding to the SU(A’c) Yang-Mills action 
the term (in Minkowski space). 
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Here ^ is a (non-negative) GF parameter, are the so- 
called Nakanishi-Lautrup multipliers and A™ = A™ [A] 
is the GF condition. For all fields we write 4) = 
where f"* are the SU(A/;) generators. Going on-shell with 
the b fields, one obtains the condition ^5™ = A™ and the 




corresponding GF action, 
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^oP = ^y'd4x(A-)2. (2) 

gauges are obtained when the linear condition A™ = 
is chosen. In this case the (non-perturbative) 
gluon propagator can be decomposed according to 

Afj,^{q) = /Q )Ax(( 7 ) + (9/i9i//9 )AL(g' ). (3) 

Slavnov-Taylor identities ensure that ^ to all 

orders and, therefore, all the dynamical information is 
carried by the transverse form factor A^. 

The lattice formulation of Yang-Mills theories is ob¬ 
tained in terms of the Wilson gauge action, in which the 
dynamical variables are the gauge links U^, related to the 
gauge fields (in lattice units) through, 

Uf^{x) = eyip[igoAf_t{x + e^/2)], 


traceless 

where is the unit vector along the direction /i, /3 = 
2Nc/gQ is the lattice coupling which determines the lat¬ 
tice spacing a. Physical quantities are then obtained by 
the evaluation of the Euclidean path integral through 
Monte Carlo techniques, with a probability distribution 
given by the exponential of the action. 

In gauges, besides the usual integration over the 
link variables, one has to integrate over the A fields. 
Eq. ([^ implies that the integration measure is a Gaussian 
distribution, with variance 

P[A™(x)](xexp|-l^ [A”^(^)]"|. (5) 

The numerical difficulty of implementing the gauges 
lies in enforcing the GF condition. 

In fact, the standard procedure for GF requires to 
gauge rotate all link variables through the gauge trans¬ 
formation U^{x) — )■ g{x)Ufj_{x)g^{x + e^), where g are 
elements of the SU(Ac) gauge group that minimizes a 
suitable functional implementing the desired GF condi¬ 
tion. In the Landau gauge case, which is the ^ 0 limit 

of the R^ gauges studied here, the functional is, 

£i.G[U,g] = -ReTr '^g{x)U^{x)g^{x+ e^), (6) 

X,fl 

which directly leads to the condition V-A’" = 0. Con¬ 
trary to this simple limit, the general case of a non¬ 
vanishing ^ was proven to have no suitable GF functional 
to minimize |44| . This no-go theorem has been evaded 
in m, where it was shown that the functional, 

^R^P.g] = Si^Gp^g] -kReTr '^ig{x)A{x ), (7) 


A^(x -k e^/2) = 


UrP) - uli 



Figure 1: (Top) The 32^ values of V-A^ evaluated for a 
configuration gauge fixed at ^ = 0.5, grouped in 5000 
bins, compared with a Gaussian with standard 
deviation ~ 0.316. (Bottom) Plot of d = V-A"^ — A'^; 
the two distributions coincide within '/d precision. 


yields the correct condition V • A"* = A™. In prac¬ 
tice, the gauge transformation g is built as a prod¬ 
uct of a sequence of infinitesimal gauge transformations 
9 = T\j^9j- For each infinitesimal transformation 6gj 
one minimizes the functional 0; however when moving 
on to the next infinitesimal transformation the 

Gaussian distribution A™ is maintained unchanged and 
the link is updated through a gauge rotation. Writing 
5gj = 1 -k J 1 variation of the functional Q 

with respect to the coefficients w"* reads. 


[t/, 5g] - [t/, 1] = Tr ^ w”^{x) ^ A(a:) , (8) 

X, m 

= ^9o A^{x -k e^/2) - Af,{x - e^/2) - A{x) . 


Thus, choosing w™ = aA™, with a is a relaxation pa¬ 
rameter to be optimized, will reduce A. Our goal is to 
converge to a vanishing A in all lattice points x, which 
implies. 


X 

Whenever this condition is fulfilled (which, based on the 
experience of Landau GF jM], means to have 9 < 10“^® 
), then the configuration is R^ gauge fixed. 

Lattice setup and algorithm. In order to study the 
gluon propagator we use 50 configurations generated 
through importance sampling of the SU(3) Wilson action 
m- We opt for a symmetric lattice of size L = 32 and 
(3 = 6.0, with associated lattice spacing a = 0.1016(25) 
fm measured from the string tension |56| . The simulated 
volume is therefore (3.25 fm)^, large enough to resolve 
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Figure 2: (Left) The longitudinal dressing function (/^Al = ^; a fit of the data to a constant yield ^ = 0.103(2), 
0.203(2), 0.302(3), 0.402(3) and 0.502(3) respectively. (Right) The R^ gluon transverse dressing function 
Landau gauge results obtained for a symmetric lattice of L = 80 and (3 = 6.0 (gray crosses) are also plotted |17| . 


the onset of nonperturbative effects in the propagator’s 
transverse form factor. The values of ^ chosen are ^ = 
0.1, 0.2, 0.3, 0.4, and 0.5. For comparison, we also report 
the Landau gauge, equivalent to Rq. 

To gauge fix the configurations, we generate, for ev¬ 
ery lattice site x, 8 real valued functions A™ using the 
Gaussian probability distribution for the generation 
of the Gaussian distribution we used the standard Box- 
Muller algorithm. After combining these functions in the 
SU(3) algebra element A, one proceeds to minimize the 
GF functional ([^. Notice that for each gauge configu¬ 
ration, we integrate over 50 different A’s. To reduce the 
correlations in the evaluation of the path integral, the 
A’s are generated independently for each gauge configu¬ 
ration. 

For minimization purposes we first tried to apply three 
different standard optimized techniques used in the Lan¬ 
dau case [^: the Fast Fourier Transform - accelerated 
steepest descent (FFT), Over Relaxation (OVR) and 
Stochastic Relaxation (STR). Each one of these three 
techniques showed a typical GF success rate of around ~ 
75% for ^ = 0.3, with the convergence rate dropping to ^ 
40% for ^ = 0.5. Thus we opted to cycle through all con¬ 
vergence techniques when the procedure stalls. Indeed, 
by cycling through FFT, OVR and STR, for our hard¬ 
est case of ^ = 0.5, we increase the convergence success 
rate up to ~90%; for the remaining 10% cases, restarting 
the combined algorithm, after performing finite random 
gauge transformations, leads to convergence for all cases 
EH- Fig-HI compares the distribution obtained from the 
values of V-A™ with the one expected for A"* for a given 
configuration and a given color index (m = 4); as one can 
appreciate the GF is within the precision defined above. 

Simulation results. The lattice gluon two-point corre¬ 
lation function reads 

{A^iq)A:(q')) = S^^A^,(q)L^S(q + q'), (10) 


where is given in Eq. (l3|. The lattice momenta 
q (used for Fourier transfornis) and q (the ‘continuum’ 
momentum) are defined according to |171 [55] , 


2 . 9m 
9m = -sm—; 
a 2 


9m = 


2'Knfi 

IT 


= 1 , 2 ,. ..,L. ( 11 ) 


From Eqs. ([^ and ( |Io| ) it follows that the transverse 
and longitudinal SU(3) propagator form factors can be 
estimated using. 


At(9^) = ^ E (<5m^ 

AT(o) = ^E(^r(o)AM(o)), 

/i,m 

^ E 9M9./9^(7i;r(9)^"(-9)). (12) 


To begin with, we show in Fig. the gluon dressing 
functions. ’Within statistical fluctuations, the longitu¬ 
dinal q^Ai^{q^) should be a constant function, coinciding 
with the variance ^ of the A probability distribution. This 
is evidently true for all cases analyzed. In the right panel 
of Fig.j^we plot the transverse 9 ^At( 9 ^) for the different 
^ values studied, renormalized at p, = 4.317 GeV using a 
momentum subtraction scheme m- Glearly, no signifi¬ 
cant deviation from the Landau gauge case is observed. 

Next, we turn our attention to the transverse form fac¬ 
tor At (Fig. left panel). As was already the case in 
the Landau gauge, one can see that the R^ transverse 
propagators show an inflection point, implying that the 
associated spectral density is not positive definite; as this 
violates the quantum field theory axiom of reflection pos¬ 
itivity, it is a manifestation of confinement m- In addi¬ 
tion, they have a marked tendency to flatten towards the 
small momentum region, thus providing strong evidence 
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Figure 3: (Left) The transverse propagator renormalized &t fi = 
I right HZ]) provide an estimate for the volume effects expected at 


4.317 [GeV]. The gray crosses (same as Fig. 
= 0. {Right) The ratio AT((i'^)/A|“°(( 7 ^). 



Figure 4: The reconstructed dynamically generated 
gluon mass in the R^ gauges. 


that also in the ^ 7 ^ 0 case the behavior of the zero- 
momentum modes of the lattice gluon field are tamed by 
the dynamical generation of a (momentum-dependent) 
gluon mass. The data confirms an IR hierarchy such that 
At (slightly) decreases for increasing values of the gauge 
fixing parameter |481150] . This is better seen in Fig. 
(right) where we plot the ratio of the transverse propa¬ 
gator to the Landau gauge propagator A|“° as a function 
of the momentum for the two values ^ = 0.1 and ^ = 0.5; 
observing a maximum difference of about 10 %. 

Comparison with continuum studies. Fig. [2 (right) 
turns out to differ from the effects observed in [ 0 I](*.e.,a 
dressing function in which the height of the peak rapidly 
increases and its location moves towards higher q^ values 
with increasing ^). Thus the effects reported are proba¬ 
bly artefacts of the truncations employed. On the other 
hand, our data agrees qualitatively with the Nielsen iden¬ 
tities analysis of Ref. (55] , where the results for the gluon 


propagator were interpreted in terms of the presence of 
a dynamically generated gluon mass. 

We now briefly apply the analysis of Ref [62] to our R^ 
lattice propagators. First we write [51] A”^ = q^J^+w?, 
where, for small one has, 


= 


“(0 + c(0 ( log ^ 


mUiq^). (13) 


To lowest order a(^) = 1-1- oi^, c(^) = Cni^, and is 

the Landau gauge dynamical gluon mass. By repeating 
the analysis of m using as input our Landau gauge L = 
32, /3 = 6.0 data, we estimate Cni ~ 0.32 [63]. Next, 
we simply fit ai by requiring that the resumed mass (131 
equals the value of At^(O) for the corresponding value 
of the GF parameter ^ and solving the renormalization 
group improved gluon mass equation [3] (yielding 
We obtain oi « 0.26, and the resulting mass is plotted 
in Fig. [^ 

Conclusions and outlook. In this paper the lattice 
SU(3) gluon propagators in R^ gauges is computed for 
large ^’s and for a lattice volume large enough to ac¬ 
cess the IR dynamics. From the numerical point of view, 
the most intensive task was the gauge fixing due to the 
large number of GF’s required and algorithmic issues. At 
least for the set of parameters simulated here, a proper 
combination of various methods solved the minimisation 
problem associated with the GF in R^ gauges; this will 
be presented in detail elsewhere HZ]. 

Our propagators show very similar characteristics 
to the one in the Landau gauge: an inflection point in 
the few hundreds MeV region followed by saturation to a 
finite value in the IR. Comparing with very recent contin¬ 
uum analytic studies, our propagators are in agreement 
with Ref. [62], which allowed us to estimate the dynam¬ 
ically generated R^ mass. This analysis suggests that 
dynamical gluon mass generation is a common feature of 
all R^ gauges in SU(3) Yang-Mills theories. 
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